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ABSTRACT 

Recent X-ray observations of galaxy clusters suggest that cluster populations are bimodally dis- 
tributed according to central gas entropy and are separated into two distinct classes: cool core (CC) 
and non-cool core (NCC) clusters. While it is widely accepted that AGN feedback plays a key role in 
offsetting radiative losses and maintaining many clusters in the CC state, the origin of NCC clusters 
is much less clear. At the same time, a handful of extremely powerful AGN outbursts have recently 
been detected in clusters, with a total energy ^ 10^^ — 10^^ erg. Using two dimensional hydrodynamic 
simulations, we show that if a large fraction of this energy is deposited near the centers of CC clusters, 
which is likely common due to dense cores, these AGN outbursts can completely remove CCs, trans- 
forming them to NCC clusters. Our model also has interesting implications for cluster abundance 
profiles, which usually show a central peak in CC systems. Our calculations indicate that during 
the CC to NCC transformation, AGN outbursts efficiently mix metals in cluster central regions, and 
may even remove central abundance peaks if they are not broad enough. For CC clusters with broad 
central abundance peaks, AGN outbursts decrease peak abundances, but can not effectively destroy 
the peaks. Our model may simultaneously explain the contradictory (possibly bimodal) results of 
abundance profiles in NCC clusters, some of which are nearly flat, while others have strong central 
peaks similar to those in CC clusters. A statistical analysis of the sizes of central abundance peaks 
and their redshift evolution may shed interesting insights on the origin of both types of NCC clusters 
and the evolution history of thermodynamics and AGN activity in clusters. 

Subject headings: cosmic rays - galaxies: clusters: general - galaxies: active - X-rays: galaxies: 
clusters - galaxies: abundances 
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1. INTRODUCTION 

The hot gas in galaxy clusters emits prolifically in X- 
rays and has been extensively studied by X-ray telescopes 
Chandra and XMM-Newton, which reveal a striking bi- 
modality in the properties of cluster cores. According 
to their core gas entropy, clusters are separated into two 
distinct classes: cool core (CC) and non-cool core (NCC) 
clusters. While the former have a low central gas entropy 

peaked at S' = k^T/v^^'^ ~ ISkeVcm^, the latter usually 
have high-entropy cores peaked at S* ^ 150keVcm^, and 
there is a disti nct less-populated gap between ~ 30 and 
~ SOkeVcm^ (jCavagnolo et al.l [20091 ). This bimodality 
also appears in cluster temperature profiles, which de- 
crease significantly toward the center in CC clusters, but 
rema in relatively flat in NCC clusters ([Sanderson et al.l 
I2006D . Observations also suggest that AGN activity 
and star formation in cluster central domi nant galaxies 
are much more pronounced in CC clusters ([Burns! [19901 : 
iCavagnolo eFa l. 2008; Raffe rtv et aT1[200l . 

The origin of this dichotomy has been studied by many 
authors. Galaxy clusters may naturally reach the CC 
state due to the interplay between radiative cooling and 
heating by active galactic nuclei (AGNs). It is widely 
thought that AGN outbursts play a key role in heating 
the intracluster mediu m (ICM), thus preventing cooling 
catastrophe (see iMcNamara k, Nu lscn 2007i for a recent 
review). It was also shown that AGN may operate as a 

^ UCO/Lick Observatory, Department of Astronomy and Astro- 
physics, University of California, Santa Cruz, CA 95064, USA; 
fulai@ucolick.org 



self-regulating feedback mechanism, which is essential in 
suppressing global thermal instabil ity and thus in m ain- 
taining the ICM in the CC state ([Guo et al.ll2008( ). In 
contrast, the origin of NCC clusters is much less clear; 
three competing explanations have been proposed: 

1. Mergers: Mergers arc shown to significantly disturb 
cluster CCs and mix the ICM (Ri cker fc Sarazp 
[2001UG6mez et al.ll2002l: [Ritchie fc Thomas[[200ir 
However, recent simulations bv [Poole et al.l ([2008[ ) 
find that CC systems are remarkably robust and 
can be disrupted only in direct head-on or mul- 
tiple collisions; even so, the resulting warm core 
state is o nly tr a nsien t. Cosmological simulations 
by B urns et al.l ([2008D suggest that NCC clusters 
may form when they experience major mergers 
early in their evolution which destroy embryonic 
CCs. Note that clusters in these simulations suf- 
fer from the overcooling problem since they do not 
incorporate mechanisms (such as AGN feedback) 
to stop cooling catastrophe. Furthermore, the rel- 
atively low numerical resolution in these cosmolog- 
ical simulations (15.6 /i^^kpc) may preclude firm 
conclusions about core structure and evolution. 



2. Pre-heating: IMcCarthv et al.l ([20081 ) suggested 
that early pre-heating prior to cluster collapse 
could explain the formation of CC/NCC systems, 
which receive lower/higher levels of pre-heating. A 
possible concern in such scenarios is whether one 
can pre-heat the ICM to a high adiabat and yet re- 
tain sufficient low entropy gas in lower mass halos 
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to obtain a realistic g alaxy population. Further - 
more, observations by iRossetti fc Molendil (|20100 
find that most NCC clusters host regions with low- 
entropy, but relatively high metallicity gas, which 
are probably a signature of recent CC to NCC 
transformation, apparently inconsistent with both 
the scenarios proposed bv [McCarthy et aTl ()2008f ) 
and lBurns et al 
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3. AGN outbursts: Strong AGN outbursts with 
^agn ^ 10^^-10^^ erg have been de tected in clus- 
ters, e. g.. Hydra A (iNulsen et al. 2005a), Hcr- 
cules AlNulsen et al.l[2 005bf) and MS0735.6-^742 1 



(iMcNamara et al.l 1200 5: McN amara et al.l l2009f ). 
Recently in lGuo fc OhI (|200a hereafter GO09), we 
show that such strong AGN outbursts could bring 
a CC cluster to the NCC state, which can be stably 
maintained by conductive heating from the cluster 
outskirts. AGN outbursts may also drag magnetic 
field lines radially, thus enhancing thermal conduc- 
tion, which could significantly heat CCs in high- 
temeprature clusters. 

The origin of NCC clusters may not necessarily explain 
the origin of the CC/NCC bimodality. To ensure the pro- 
nounced bimodality seen in observations, clusters need to 
stay in the NCC state for a duration at least comparable 
to the cooling time. In fact, many NCC clus ters have a 
short cooling time (~ 1 Gvr: [Sanderson et al.i r2006'). and 
the NCC peak in the cluster central entropy distribu- 
tion may not exist unless heating largely offsets radiative 
cooling for a relatively long time. While clusters may be 
maintained in the CC state by episodic AGN feedback, 
thermal conduction can stably keep clusters in the NCC 
state (Guo et al. 2008). This has been confirmed by nu- 
merical simulations in GO09, which further predicts that 
the CC/NCC dichotomy is more pronounced in higher- 
temperature clusters, due to the fact that the heating 
and stablizing effects of conduction decline with temper- 
ature. 

In GO09, we adopted a simplified 'effervescent model' 
for AGN heating (Begelman 2001) . and performed one- 
dimensional (ID) calculations. X-ray cavities in this 
model are in pressure equilibrium with and rise buoy- 
antly in the ICM, resulting in a very gentle AGN heat- 
ing. Furthermore, the ID model assumes spherical sym- 
metry. However, real AGN outbursts produce jets and 
X-ray cavities in two opposite directions, which are by no 
means spherically symmetric. The formation and evolu- 
tion of X-ray cavities are probably much more dynamic 
than assumed in the 'effervescent model'. Shock waves 
have been detected in many clusters with X-ray cavities, 
and are usually thought to be a natural result of AGN en- 
ergy released in the ICM. In this paper, we study the evo- 
lution of X-ray cavities formed as AGN energy is injected 
into the ICM, using two-dimensional (2D) hydrodynam- 
ical simulations. Our primary goal is to investigate if 
strong AGN outbursts can transform a CC cluster to the 
NCC state, and if so, how this transformation happens. 

The metallicity distribution of the ICM, extensively 
measured with X-ray observations, contains important 
clues about the physics of galaxy clusters. The mass 
and distribution of metals in the ICM constrain the inte- 
grated history of past star formation (metals are released 



via supernovae explosio ns and winds) and the ICM en- 
richment processes (see iBohringer fc Werner! 120101 for a 
recent review). The spatial abundance distribution is 
also significantly affected by transport processes in the 
ICM, e.g., turbule nt mixi ng triggered b y central AGN 
outbu rsts (R ebusco et al.l [2005.. 20061 iRoediger et al.l 
120071 ). Of great interest is the bimodality in cen- 
tral abundance profiles between CC and NCC clus- 
ters. CC clusters usually have a peak in the iron abun- 
dance profile at the cluster center, while many NCC 
clus ters show a relatively flat radial abundance pro- 
file (IDe Grandi fc Molendil IMlt IDe Grandi et al.l[200l 
iLeccardi et al.l 20101). Neverthel ess, recent observa- 
tions b y lLeccardi fc Molendl ()2008D and lSanderson et all 
(|2009f ) suggest that some NCC clusters also have central 
iron abundance peaks. The differences in abundance pro- 
files between CC and NCC systems have not been stud- 
ied with detailed computational models. If NCC clusters 
are transformed from CC systems (i.e., not primordial), 
the same process transforming the CC to NCC systems 
may also explain the differences seen in their abundance 
profiles. In this paper we follow the evolution of iron 
abundance in our simulations and investigate if strong 
AGN outbursts can significantly change abundance pro- 
files as they transform a CC cluster to the NCC state. 
We show that strong AGN outbursts efficiently mix met- 
als in cluster central regions, but the final abundance 
profile in the NCC state strongly depends on the size 
of the initial abundance peak at the CC cluster cen- 
ter. Our model may naturally explain the the range of 
abundance profiles observed in NCC clusters. CC sys- 
tems with broad/narrow central abundance peaks can 
be transformed by powerful AGN outbursts into NCC 
systems with/without central abundance peaks. A sta- 
tistical analysis of the sizes of abundance peaks and their 
redshift evolution could test the validity and applicabil- 
ity of our model and may shed interesting insights on the 
origin of the CC/NCC bimodality. 

The rest of the paper is organized as follows. In Sec- 
tion [21 we present our model, including basic equations 
and numerical setup. Our results are presented in Sec- 
tion |3l We summarize our main results in Section |4| with 
a discussion of the implications. The cosmological pa- 
rameters used throughout this paper are: rim = 0.3, 
JIa = 0.7, h — 0.7. We have rescaled observational re- 
sults if the original paper used a different cosmology. 



2. EQUATIONS AND NUMERICAL METHODS 

2.1. Basic Equations 

AGN outbursts inject cosmic rays (CRs) into the ICM, 
producing X-ray cavities in the cluster gas which have 
also been observed at radio frequencies in many clusters 
due to synchrotron emission of relativistic electrons. The 
outbursts induce weak shocks, and drive l arge amounts of 
gas mass outflow HGuo fc Mathew!j|2010[ ). In this paper, 
we study the combined evolution of thermal gas, CRs, 
and the gas metallicity during and after outbursts. The 
governing equations may be written as: 



dp „ 



0, 



(1) 
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P^ = -V(F + Pc)-pV$, 



(2) 



3f 

— + V-(ev) = -PV-v-nineA(r,Z), (3) 
ot 



dec 
dt 



V-(ecV) = -PeV-V + V-(KVee) + S'c, (4) 



dt 



- PFcV ■ 



(5) 



where d/dt = d/dt + v • V is the Lagrangian time deriva- 
tive, Pc is the CR pressure, Cc is the CR energy density, 
K is the CR diffusion coefficient, Sc is the CR source term 
due to the central AGN activity, pp^ is the iron density, 
and all other variables have their usual meanings. Pres- 
sures and energy densities are related via P = (7 — l)e 
and Pc = (7c — l)ec, where we assume 7 — 5/3 and 

Equation ([5]) describes the conservation of iron mass. 
Since we focus on the effect of AGN outbursts on the 
iron distribution and follow the cluster evolution for a 
timescale much shorte r than enrichment times (^5 Gyr; 
lB5hringer et al.ll2004D . we ignore the iron source term. 
The iron abundance Z in units of the solar value is pro- 
portional to Z cx pFe/p- Thus the iron density ppo in 
equation ((5]) may be replaced by Zp. Since both the gas 
mass and iron mass are conserved, the metallicity Z is 
also conserved: 

dZ__dZ_ 
H ^ ~dt 



V • VZ = 0, 



(6) 



and therefore provides a tracer for the ICM gas, which is 
helpful in understanding how AGN outbursts affect and 
mix the ICM. 

In the gas energy equation ([3]), we include radia- 
tive cooling with a volume cooling rate nincA{T, Z), 
where the cooling fu nction A(r, Z) is adopted from 
[Sutherland fc Dopital ([1993) and depends on both gas 
temperature T and metallicity Z . The ion number den- 
sity Hi is related to the proton number density nn via 
rij = l.lnu, and thus the molecular weight is p = 0.61. 
The gas temperature is related to the gas pressure and 
density via the ideal gas law: 



T = 



kBP 



(7) 



where fee is Boltzmann's constant and is the atomic 
mass unit. 

Equations ([l]) — ([5]) are solved in (r, z) cylindrical co- 
ordinates usin g a two-dimensional Eu lerian code similar 
to ZEUS 2D (jStone fc Normanlll99^ : in particular, we 
have incorporated into the code a background gravita- 
tional potential, CR diffusion, CR energy equation, and 
iron equation of mass conservation. The computational 
grid consists of 100 equally spaced zones in both coordi- 
nates out to 100 kpc plus additional 100 logarithmically- 
spaced zones out to 1 Mpc. For all the three fluids, we 
adopt reflective boundary conditions at the origin and 
outflow boundary conditions at the outer boundary. 
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Fig. 1. — Initial radial profiles of electron number density (top 
panel) and temperature (bottom panel). The dotted line in tfie 
bottom panel is th e be st analytical fit to Chandra data between 
40 to 1000 kpc by IVlkhlinin et al. (2001 ); cr osses in both panels 
indicate Chandra data bv iEttori et alTl2002l) . Our initial cluster 
profiles (solid lines) provides a very good fit to data for the whole 
cluster out to 1 Mpc. 
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Fig. 2. — Initial abundance profiles. The solid line (adopted in 
our main run Dl) has a central abundance peak with size Tnictal = 
160 kpc and fits Chanda data of A1795 (crosses; Ettori ct al. 2003) 
reasonably well, which have large scatter. The dotted line has a 
much smaller central abundance peak (r^ietal = 60 kpc) and is 
adopted in ru n D2. Chandra data for the cluster A2199 (triangles; 
I Johnstone et al. 2002), which has a much smaller central abun- 
dance peak, is shown for comparison. 



2.2. Initial Cluster Profiles 

Our model and methods are generally applicable to all 
relaxed clusters, but for concreteness, we adopt simula- 
tion parameters appropriate for the typical CC cluster 
Abell 1795, which has b een well observed by both Chan- 
dra and XMM-Newton (iTamura et al.|[200ltlEttori et al.l 
[200l lVikhlinin et eiU\20Ml 

For initial profiles of A1795, we ffist build an analytic 
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fit to the deprojected 3-dimensional gas temperature pro- 
file derived from Chandra observations, which acquired 
data out to ^ 1 Mpc covering our entire computational 
grid: 



T(r) = (TiS + r^) 



l/a 



(8) 



where the constant a is chosen to be a = 5, Tin = 2.5 
keV is the observed central temperatur e of A1795, and 
Ty is the best- fit temperature profile of iVikhlinin et afl 
((2006) which provides an excellent fit to Chandra data 
of A1795 from 40 kpc to 1 Mpc: 

(9) 



where 



rv(r) =9.68tcooi(r)tout(r) keV, 



tr 



lir) 



1.03 



0.1 + (r/ri) 
1 + (r/ri)i-03 



(10) 



describes that the temperature declines inward in the 
central region, and 



tout{r) 



1 



[1 (r/ra)"]''/'^ ' 



(11) 



represents the radially declining outer region. Here the 
parameters are ri = 77 kpc, r2 = 550 kpc, a = 1.63 and 
b = 0.9. The resulting temperature profile T{r) is shown 
as the solid line in the bottom panel of Figure [H where 
T\f(r) (the dotted line) and Chandra data (crosses) from 
lEttori et al.l ()2002[ ) are also shown. Clearly our analytic 
temperature profile provides a very good fit to observa- 
tions out to ~ 1 Mpc. 

The initial cluster density profile is solved by assuming 
hydrostatic equilibrium. At the beginning of our simu- 
lation t = 0, the CR energy density is assumed to be 
Ccr — throughout the cluster. The gravitational poten- 
tial $ is determined b y the dark matter N avarro-Frenk- 
White (NEW) profile ([Navarro et al.llT997[) : 



PBuir) 



Mo/2'K 
r{r + rsY ' 



(12) 



where r^ is the standard scale radius of the NEW profile 
and Mq is a characteristic mass. This density distribu- 
tion results in a gravitational potential: 



$ = - 



2GMoln(l-Kr/rs) 

Ts r/ rs 



(13) 



Eor A1795, we choose Mq = 4.8 x lO^'^Mp, an d r^ = 
460 kpc. At redshift 0.0632 ()Ettori et al.l 1200^ . these 
parameters correspond to a virial mass Mvir = 7.87 x 
IO^^'Mq, and a concentration c — 4.07 assuming that the 
virial radius is Tvir r2m^ within which the mean mass 
density is 200 times the critical density of the Universe. 
Taking the electron number density rig at the origin to be 
no = 0.065 cm~3, we derive from hydrostatic equilibrium 
the initial radial profile of ric, which fits observations very 
well, as clearly shown in the top panel of Figure [1] 

X-ray observations have been extensively used to 
measure abundance profiles in galaxy groups and 
clusters. In CC clusters, the observed abundance 
distributions typically peak in the central, low- 
entropy core regions, and d ecline radially o u tward 
(|De Grandi fc Molendil [200lL IDe Grandi etall [200l 



iBaldi et al.|[2007t iLeccardi fc Molendill2008[ ). These ob- 
servations also suggest that abundance profiles flatten 
off with Z ~ 0.2 — 0.3 at large radii. In this paper the 
metallicity Z is calculated relative to the solar iron abun- 
da nce (ZfT) =Ee/H= 4 . 68 x 10~^ by number) published 
bv lAnders fc Grevessd (jl989() . It has been superseded by 
the ne w va lue Zp-, = 3.16 x 10~ ^ of IGrevesse fc Sauvall 
(|1998[) and lAsplund et al.l (|2009[ ). but it allows straight- 
forward comparison with most of t he literature. Further- 
more , the metallicity Z used in the iSuther land fc Dopital 
(11993) cooling func t ion A jT, Z) is also in units of the 
lAnders fc Grevessd ()1989[ ) iron abundance. A simple 
scaling by 1.48 converts our abundance values to the m et- 
alicities relative to the new IGrevesse fc Sauvall (|1998l ) so- 
lar abundance, without changing any other result in our 
paper. 

The initial abundance profile is chosen to be 



z = {z'^ + z^Y 



in 



(14) 



where /3 = 5. The constant Zq = 0.27 sets a minimum 
for the abundance profile, which repre sents the fiatten- 
ing off at large radii, and is taken from IDe Grandi et al.l 
(i20(]4|'). wh ere we have converted metalicities in unit of 
the IGreves se & Sauval ( 19 9j) iron a bunda nce to those 
relative to the lAnders fc Grevessd (1989) abundance. 
The outward abundance decline is described by 



= 0.1 



(15) 



where rmctai characterizes the spatial size of the central 
metallicity peak. The initial abundance profiles for our 
runs and Chandra data are shown in Figure [51 

2.3. Cosmic-ray Physics and Assumptions 

X-ray cavities are usually thought to be infiated by 
bipolar jets emanating from AGNs located at cluster 
centers. Radio synchrotron emission has been detected 
from many X-ray cavities (especially when the cavi- 
ties are young), confirming the existence of a signifi- 
cant amount of CRs, which may be transported along 
the jets and/or created in strong shocks as the jets en- 
counter the ICM. One can envision that the jets de- 
posit CRs into small regions at their terminal points, 
which expand and form underdense bubbles producing 
the observed X-ray cavities. We have thus developed 
a model to numerically study CR injection and evo- 
lution in the I CM, which has been successfully used 
in Mathews fc BriEhcnti (2008a), Mathews fc Brighcirt^ 
(|2008bD . IMathewsl '(2009). and Guo fc Mathewsl l\20m . 
where the reader is referred to for further details. Here 
we simply summarize several modifications and reiterate 
some important points. 

The injection of CRs into the ICM is described in equa- 
tion ^ by the source term Sc- We assume that the CRs 
are deposited into a Gaussian-shaped sphere of charac- 
teristic radius rg = 5 kpc located at rcav = (^i -z) = 

(0, ^cav): 



^3/2^3 







when t < iagn, 
when t > iagn, 



(16) 



where tagn = 10 Myr is the duration of the CR injection 
(AGN active phase), and i^agn is the total injected CR 
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TABLE 1 
List of Simulations 



-C/agn ^cav ' metal ' 

Run (10*^^ erg) (kpc) (kpc) (kpc) 



Dl 


3.15 


10 


160 




Dl-A 


3.15 


10-50'= 


160 




Dl-B 


6.3 


10 


160 




Dl-C 


3.15 


30 


160 




Dl-D 


3.15 


10 


160 


10 


D2 


3.15 


10 


60 




D2-A 


3.15 


10 


40 





"The AGN energy released in the form of CRs in one 
hemisphere during the active phase t < tagn = 10 Myr. 

'i"cav = (0, ZcB.v) is the position where the cosmic rays 
are injected during the AGN active phase. 

'^''motal characterizes the sp atial size of the central 
abundance peak (see equation 11511 . 

'^ro is an inner cutoff radius characterizing a central 
dip in the initial cluster abundance prof ile, w hich is 
only introduced in run Dl-D (see Section l3.2| l. 

■^In this run, Zcav moves from 10 to 50 kpc with a 
constant speed during the CR injection phase. 

energy in one bubble (2i?agn for the whole cluster). The 
integral of over space gives the the CR injection lu- 
minosity i?agn/iagn during the active phase. See Tabled] 
for specific model parameters in each run. 

In addition to their advection with the thermal gas, 
CRs diffuse through the gas as described in equation (|4]). 
The CR diffusion coefficient k is poorly known but may 
vary inversely with the gas density sin ce the magneti c 
field is probably larger in denser gas ((Dolag et al.ll2Q0l[ ). 
We adopt the following functional dependence of the dif- 
fusion coefficient on the gas density: 

_ riO^°(noo/no) cm^ s"^ when no > rieo, 
\lO'^° cm^ s~^ when no < noo, 

where noo — 10^^ cm^"^. This level of diffusion does 
not strongly affect the early evolution of cavities, and 
our results are fairly insensitive to it. During their dif- 
fusion, CRs interact with magnetic irregularities and 
Alfven waves, exerting CR pressure gradients on the 
thermal gas (equation [2]). The early evolution of cav- 
ities and associated shock heating are mainly deter- 
mined by the expansion of the surrounding gas driven 
by the CR pressure. We neglect other more complicated 
(probably secondary) interactions of CRs with thermal 
gas, e.g.. Coulomb interactions, hadronic collisions, and 
hydromagnetic-wave-mediated CR heating, that depend 
on the CR energy spectrum and provide addition al heat- 
ing effects for the ICM (e.g., IGuo fc OhI 12001 . Our 
model reproduces very well many characteristic morpho- 
logical features of the huge X-ray cavities observe d in the 
cluster MS 0735.6+7421 (|Guo fc Mathewsl [20T0I ) which 
experienced an AGN energy release comparable to that 
considered here. 

3. RESULTS 

3.1. Removal of Cool Cores 

The primary goal of this paper is to determine if pow- 
erful AGN outbursts can transform CC clusters to the 
NCC state, and if they can, how central abundance peaks 
typically found in CC clusters are affected during this 




r (kpc) 

Fig. 3. — Time evolution of radial profiles of emission- weighted 
spherically averaged (a) electron number density, (b) temperature, 
and (c) entropy in our main run Dl. After t ~ 0.2 — 0.3 Gyr, the 
cluster relaxes to the NCC state with a relatively flat temperature 
profile and a high central entropy (~ 100 keV cm^) core. 

transformation. We begin with the removal of cool cores 
in this subsection. Our model is intended to be generic, 
but for definiteness, we choose the CC cluster A1795 as 
our fiducial cluster. Model parameters in our runs are 
listed in Table [TJ We first present the results of our main 
run Dl, and then compare this with runs computed with 
other combinations of the parameters in Table 1. 

Figure [3] shows time evolution of radial profiles of 
emission-weighted spherically averaged gas properties, 
including (a) electron number density, (b) temperature, 
and (c) entropy in run Dl. The solid lines represent the 
observed profiles of the CC cluster A1795, which serve 
as the initial state in our simulations. Within the CC 
(< 100 kpc), both the temperature and entropy decrease 
inward and reach low values T ^ 2.5 keV, and S ^ 10-20 
kev cm^ at the cluster center. As AGN energy in the 
form of CRs is injected into the ICM, shocks appear and 
heat the gas, as clearly seen in the top-left panel of Fig- 
ure m (the shock front is located near the outer edge of 
the red 'annulus'), which shows central slices (100 x 100 
kpc) of electron number density, entropy, and metallicy 
distributions at three times t=10 (top panels), 100 (mid- 
dle panels), 200 Myr (bottom panels). The dotted line 
in Figure [Sl^c) suggests that all the gas in the CC has al- 
ready been heated to a high entropy ^ 100-200 kev cm^ 
at t = 50 Myr, but Figure [3ja) also indicates that at this 
time CRs have created a huge low-density cavity at the 
center. 

The formation and final break-up of the cavity may be 
best seen in Figured] The top and middle panels show 
that the cavity is initially formed along the z direction 
where CRs are injected, and then merges with the lower 
cavity (z < 0; not shown in Figure [5]) at the reflective 
boundary z — 0. The cavity is formed as the CR pres- 
sure causes the expansion of the surrounding gas, driv- 
ing a large gas outflow, as indicated in Figure [S] which 
shows the time evolution of the gas mass Mgas(< r) en- 
closed within various radii r = 10, 20, 50, 100 kpc (see 
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Fig. 4. — Central slices (100 X 100 kpc) of log (?ic/cm~'^) (left column), log (S/keV cm^) (middle column), and log (Z/^solar) (right 
column) at three time epochs t=10 (top), 100 (middle), 200 Myr (bottom) in run Dl. The horizontal and vertical axes represent r and 
z in kpc respectively. Cavity regions in the left and middle panels are saturated to show the variations of gas density and entropy in the 
ICM. The images clearly show the formation and final break-up of the low-density cavity produced by the AGN outburst. 

line in Figure [DJc) . We note that the amount of shock- 
heated gas within the cavity is very small and negligible 
compared to the inflowing thermal gas. 

The bottom panels of Figure 2] also suggest that, as 
thermal gas flows back to the cluster center, the huge 
X-ray cavity is disrupted into two pairs of cavities in the 
orthogonal direction. This morphology is quite different 
from that of regular AGN feedback events, which usually 
produce one pair of X-ray cavities in opposite directions 
and are essential to maintain many clusters in the CC 
state. However in this paper, we investigate much more 
powerful AGN outbursts, which are much rarer and may 
transform CC clusters to the NCC state. Depending on 
where the CRs are deposited, these powerful AGN out- 
bursts may produce X-ray cavities with different mor- 



IGuo fc Mathewi[2010l for more discussions and implica- 
tions). The outflow begins at small radii, and gradually 
affects gas at larger radii. The size (radius) of the cavity 
is very large, extending out to 80-100 kpc a.t t = 100 
Myr. At later times, the cavity breaks up and thermal 
gas flows back to the cluster center as seen in the bot- 
tom panels of Figure H) Figure [5] suggests that the gas 
inflow happens during t ^ 100-200 Myr. The middle 
panels of Figure 2] show a feature of gas inflow near the 
image centers, indicating that the gas indeed starts to 
flow inward as early as t ~ 100 Myr, and the gas inflow 
is not spherically symmetric. The inflow of thermal gas 
also produces shocks, which propagate within the cav- 
ity and heat some low-density gas to very high entropy 
(a few hundred keV cm^), as shown by the short-dashed 
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Fig. 5. — Time evolution of gas mass Afgas(< r) enclosed within 
various radii r = 10,20,50,100 kpc in run Dl. The AGN out- 
burst deposits CRs into the ICM, which cause a large gas out- 
flow, creating a low-density cavity at the cluster central region. At 
t ~ 0.1 — 0.2 Gyr, the cavity starts to break up as thermal gas falls 
back to the center. 

phologies. In run Dl-A, where the CRs are injected into 
the ICM continuously by a source moving out from the 
central AGN, the resulting central cavity is elongated in 
the jet direction and finally breaks up into three pairs of 
cavities in the same direction (see Figure [7]) . Interest- 
ingly, the morphology seen in this run may have already 
been detected in cluster observations (e.g., in Hydra A; 
see Figure [7] and the corresponding discussion). 

As high-entropy thermal gas flows to the center, the 
cluster relaxes to the NCC state after t ^ 0.2 — 0.3 Gyr, 
as indicated by the long-dashed and dot-short dashed 
lines in Figure [S] During the NCC state, the cluster has 
a relatively flat temperature profile and a high entropy 
(> 80 keV cm^) core. The central cooling time in the 
NCC state is quite long (~ 3 - 4 Gyr) and the cluster 
profiles do not evolve much from t = 0.3 to 0.5 Gyr (dot- 
long dashed line). 

In our calculations, we ignore thermal conduction, 
which may remove irregularities in the temperature dis- 
tribution and may even provide a strong heating source 
for the ICM du r ing the CC to NCC transformation 
(IGuo fc Ohl(200l IRuszkowski fc Ohl[2009t IParrish"eraII 
|2010[) . Iiicluster cool cores where the gas temperature 
decreases in the direction of gravity, conduction may be 
strongly suppressed by the heat flux driven buoyancy 
instability (hereafter HBI; see Quataert 2008), which re- 
orients the magnetic field to be largely transverse to the 
radius. However, it is possible that AGN outbursts coun- 
teract the HBI by disturbing the azimuthal nature of 
magnetic fields, thus enabling conduction and lowering 
the required AGN energy to accomplish the CC to NCC 
transformation. A detailed study of this process is clearly 
beyond the scope of this paper (but see Guo fc Oh . 20091 
for a simplified one-dimensional study of this effect on 
galaxy groups and clusters). 

We conducted a parameter study of our model by per- 
forming a series of simulations with varied parameters as 
listed in Table 1. The resulting radial profiles of emission- 
weighted spherically averaged gas entropy att~0.5 Gyr 
in these runs are shown in Figure [51 In run Dl-B, the 
injected AGN energy is doubled (i?agn = 6.3 x 10^^ erg), 
resulting in much higher core entropies at t = 0.5 Gyr. 
Varying the initial abundance profile has a negligible ef- 



fect on cluster evolution, as shown in runs Dl, Dl-D, 
and D2, all of which show a very similar entropy pro- 
file a.t t = 0.5 Gyr. Obviously metallicity only affects 
the radiative cooling rate, which is not significant during 
the first one half Gyr in our simulations (shock heating 
increases the cooling time of thermal gas as well). 

However, the evolution of core entropy does depend on 
where the CRs are injected. When the CRs are injected 
at a larger radii Zcav — 30 kpc in run Dl-C, the cluster 
can not be heated to the NCC state, as clearly shown 
by the dot-short dashed line in Figure [H In this case, 
a pair of X-ray cavities are produced in opposite direc- 
tions a nd do not merge at th e cluster center as in run 
Dl Csee lGuo fc Mathewsl[20lol for more details of such a 
simulation in the cluster MS 0735.6-1-7421). Even when 
the injected CR energy is doubled, the gas within the 
central ^ 10-20 kpc still can not be heated to high en- 
tropies as in run Dl. Since this very central region in 
CC clusters is usually quite dense, it is very likely that 
AGNs often deposit quite a large fraction of energy di- 
rectly into this region. If so, powerful AGN outbursts 
can remove CCs as in run Dl. In run Dl-A, we assume 
that the CRs are injected into the ICM continuously by 
a (jet) source moving out from the central AGN (more 
specifically, Zcav increases from 10 to 50 kpc at a con- 
stant speed within tagn)- The resulting entropy profile at 
t = 0.5 Gyr (the short-dashed line) clearly indicates that 
the cluster is heated to the NCC state, and, as expected, 
has higher entropies in the region ^ 10-50 kpc than in 
run Dl. It is of great interest to investigate what deter- 
mines the position where AGN energy is deposited and 
forms X-ray cavities. Incorporating the relevant physics, 
this may be studied by high-resoluton jet simulations, 
which may also reveal some important processes that 
may affect the cavity evolution an d that are not included 
in our current sirn ulations (e.g., iSternberg et al.l 120071 : 
[Sternberg fc Soked[2008i) . 

Figure [7] shows central slices of log (uc/cm^^) at three 
epochs t = 10, 100, 200 Myr in run Dl-A. The injected 
CRs produce a low-density cavity, which is much more 
elongated in the z direction than that in run Dl. As the 
cavity breaks up, its morphology resembles three pairs of 
X-ray cavities, as seen in the right panel {t = 200 Myr). 
Such multiple pairs of AGN bubbles have been detected 
in the cluster Hy dra A (denoted a s the "Swiss-cheese- 
like" topology by I Wise et al.l I2007D , which is currently 
hosting a very powerful ^ 10^^ erg AGN outburst, prob- 
ably transforming that cluster to the NCC state. This 
also suggests that the multiple pairs of AGN bubbles in 
Hydra A are probably not an expression of the AGN duty 
cycle. 

3.2. Are Central Abundance Peaks also Removed? 

Heavy metals observed in the hot ICM using X-ray 
spectroscopy indicate that in CC clusters the iron abun- 
dance pro files have central peaks an d decline radially 
outward (De Grandi & Molendi" '2001: 'De Grandi eLaU 
2004; Baldi et al. 2007; Leccardi & Molendi 2008). Inter- 
estingly, observations reach contradictory results on the 
ab undance profiles of NCC cluste rs. It has been claimed 
bv iDe Grandi fc Molendi (|200l and iDe Grandi eTall 
(pool " that NCC clusters have a nearly uniform spa- 
tial distribution of metals. In contrast, observations by 
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Fig. 6. — Radial profiles of emission- weighted spherically aver- 
aged gas entropy at t = 0.5 Gyr in various runs. The same initial 
profile for all our runs is also shown as the solid line for comparison. 

ILeccardi fc Molendil (|2008D and ISanderson et all (|2OO90 
indicate that the abundance profiles of NCC clusters are 
very similar to those of C C clusters, both showi ng cen- 
tral peaks. More recently, ILeccardi et al.l (|2010[ ) found 
both types of NCC clusters in their cluster sample. In 
this subsection, we investigate how the centrally-peaked 
abundance profile evolve when the CC cluster is trans- 
formed to the NCC state by AGN outbursts. We con- 
struct a physical scenario that can explain both obser- 
vationally contradictory results on the NCC abundance 
profiles. 

We first investigate the temporal evolution of the abun- 
dance profile in run Dl. Figure |8] shows radial profiles 
of emission-weighted spherically averaged metallicity Z 
at various time epochs in this simulation. Clearly dur- 
ing the early times (the dotted and short-dashed lines) 
as the low-density cavity is created at the cluster center, 
the abundance profile is shifted toward larger radii, while 
the metallicity within the cavity is nearly constant. As 
the cavity breaks up and the cluster relaxes to the NCC 
state, high-metallicity gas moves back to the cluster core. 
As clearly seen in Figure[5](the dot-short dashed and dot- 
long dashed lines), the gas abundance within the central 
^ 50 kpc in the NCC state is efficiently diluted by mix- 
ing, while the abundance profile at large radii remains 
similar to that in the original CC state. However, since 
the size of abundance peak in the original CC abundance 
profile is very broad (rmetai = 160 kpc), the abundance 
profile in the final NCC state still has a central peak, 
with a maximum abundance value of around 0.6 — 0.65, 
which is roughly the abundance value of the initial CC 
profile at a cluster-centric radius of 50 kpc. Thus we find 
that in this run, though the abundance peak value does 
decrease, the central abundance peak is not removed dur- 
ing the CC to NCC transfor mation, which is consistent 
with recent observat i ons by ILeccardi fc MolendU ()2008f ) 
and ISanderson et al.l (j2009'), while not explaining NCC 
clu sters with relatively fia t abun dance profiles observed 
bv iDe Grandi fc Molendl (|200l and iDe Grandi eTaD 
(|2001k 

Figure ini^a) shows radial profiles of emission- weighted 
spherically averaged metallicity at i = 0.5 Gyr in our 
runs Dl, Dl-A, Dl-B, and Dl-C. Clearly the central 



abundance peak is not removed in any of these runs, sug- 
gesting that it is very difficult to remove it by varying the 
energy or location of the CR injection. X-ray observa- 
tions indicate that some clusters and elliptical galaxies 
exhibit a dip in abundance in their very centers, e.g., 
Abell 2199 (see the long-dashed line in Figure [2|). We 
performed another run Dl-D to study if such a CC abun- 
dance profile may result in a final NCC abundance profile 
without the central abundance peak. To implement a dip 
in abundance profile at the very center, we introduce an 
inner cutoff term to in Equation (|15p : 

Z, = 0.8e-''/'-'-""(l - e"''/'^"), (18) 

where rg = 10 kpc is a characteristic inner cutoff radius. 
The initial abundance profile in run Dl-D is then cal- 
culated using Equation (ITil) . and is shown as the solid 
line in Figure IHJb). The profile has a peak abundance 
of ~ 0.6-0.65 at a cluster-centric radius of ^ 20-30 kpc, 
which is roughly the same radius where the peak in the 
abundance profile of Abell 2199 is located. The dotted 
line in Figure |9ljb) shows the spherically averaged radial 
abundance profile at t = 0.5 Gyr when the cluster is 
in the NCC state, clearly indicating that the abundance 
peak is not removed, though the metals within ^ 50 kpc 
are efficiently mixed. This result is reasonable, since the 
total amount of gas within the central abundance dip is 
relatively small (the gas mass within the central 20 kpc 
is roughly one tenth of that within the central 50 kpc). 

We performed runs D2 and D2-A to explore the depen- 
dence of the final NCC abundance profile on the peak 
size Tniotai of the initial CC abundance profile. This 
is motivated by the fact that in some clusters r^dai 
is much smaller (e.g., rmotai ^ 40-60 kpc in A2199, as 
seen in Figure [2]) than we assumed in our previous runs. 
The results are shown in Figure [TU] For run D2 with 
fmctai — 60 kpc, the final NCC abundance profile is rel- 
atively flat and has a central value of around 0.4 in the 
cluster central regio ns, which is consis t ent with observa- 
tions by |De Grandi fc Molendil (|2001f ). IDe Grandi et al.l 
(|2004f) . and .Baldi et al.l (|2007l) . In run D2-A with a 
smaller peak size rmetai = 40 kpc, the resulting NCC 
abundance profile is even flatter with a smaller central 
value of around 0.3. 

From Figure [TOl we conclude that the final NCC 
abundance profile depends sensitively on the peak 
size Tniotai of the initial CC abundance profile. It 
can be effectively removed when rmotai is compara- 
ble to or smaller than a characteristic radius 50- 
60 kpc ^ 0.03r2oo for A1795) within which power- 
ful AGN outbursts can efficiently mix metals. In- 
terestingly, X-ray observations find two contradictory 
(possibly bimodal) types of NCC clusters according to 
their central abundance profiles: NCC clusters with 
(jLeccardi fc Molendil l2008t jSanderson et al.' 2009' ) and 
without central abund ance p eaks (Dc Grandi & M olendil 
[200ltlDe Grandi et al.l[200l . Our analysis can thus nat- 
urally explain both types of NCC clusters, which may 
be formed by powerful AGN outbursts from CC clusters 
with central abundance profiles having different spatial 
sizes. 

3.3. Implications on the Evolution History of Galaxy 

Clusters 
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Fig. 7. — Central slices (100 X 100 kpc) of log {nc/cm~^) at three time epochs t=10, 100, 200 Myr in run Dl-A. The horizontal and 
vertical axes stand for r and z in kpc respectively. The images are saturated inside the cavities to show the va riations of g as density in the 
ICM. The right image shows a "Swiss-cheese-like" topology, which has been detected in the cluster Hydra A ( Wise et aI.ll2007A . 




Fig. 8. — Radial profiles of emission- weighted spherically aver- 
aged metallicity Z at various times in our main run Dl. In the final 
NCC state, the central abundance peak is not removed, though the 
peak abundance does decrease. 

X-ray observations suggest that the iron excess asso- 
ciated with central abundance peaks is mainly produced 
by th e brightest cluster galax ies fBCGs: lDe Grandi et alJ 
120041 : iBohringer et alJ I2004D an d the enrichment time s 
are usually quite long (> 5 Gyr; IBohringer et al.|[2004[ ). 
In the absence of mixing, the metallicity profiles should 
follow the optical light profiles of the BCGs and have 
narrow central abundance peaks, easily destroyed when 
strong AGN outbursts transform CC clusters to the NCC 
state. However, recent observational and theoretical 
studies suggest that CC clusters are maintained in the 
CC state by numerous AGN feedback events, which also 
gradually broaden cent ral abundance peaks by g enerat- 
ing turbulent mixing (Rebusco et al.ll2005ll2006[) or gas 
outflows ( Guo fc Mathews 2010,) . If a cluster stays in the 
CC state for a very long duration, it would build up a 
broad central abundance peak, which is very robust to 
the CC to NCC transformation. If our explanation for 
the origin of NCC clusters is correct, it implies that the 
two types of NCC clusters have different evolution histo- 
ries. NCC clusters without central abundance peaks were 
never in a CC state for a very long duration, while those 
with central abundance peaks could have experienced a 
long CC state. 

NCC clusters may cool gradually and reach the CC 
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Fig. 9. — Radial profiles of emission- weighted spherically aver- 
aged metallicity at i = 0.5 Gyr in our runs (a) Dl, Dl-A, Dl-B, 
Dl-C, and (b) Dl-D. The solid line in the top panel represents 
the same initial profile for runs Dl, Dl-A, Dl-B, and Dl-C, while 
that in the bottom panel shows the initial abundance profile in run 
Dl-D. The central abundance peak is not removed in any of these 



state. During this transformation, the cluster gas slightly 
settles inward, which may affect metallicity profiles. In 
particular, for NCC clusters with central abundance 
peaks, the sizes of abundance peaks may decrease. We 
study this effect by following run Dl for much longer 
times. Figure [TT] shows time evolution of radial pro- 
files of emission-weighted spherically averaged gas en- 
tropy and metallicity in this run. At t ~ 2.6 Gyr, the 
cluster reaches a state very close to the initial CC state, 
and the abundance profile does slightly shift toward the 
cluster center. However, the decrease in the spatial size 
of the abundance peak is small, and rmctai (~ 130 kpc) 
at t ^ 2.6 Gyr is much larger than the characteristic 
radius 50-60 kpc) within which powerful AGN out- 
bursts can efficiently remove the metallicity peak. At 
t ^ 3 Gyr, the cluster core cools to much lower entropies, 
but the abundance profile does not change appreciably 
compared to that at i ~ 2.6 Gyr, retaining a broad cen- 
tral abundance peak. This is reasonable, since during 
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Fig. 10. — Radial profiles of emission-weighted spherically aver- 
aged metallicity at i = and t = 0.5 Gyr in our runs Dl, D2, 
D2-A. The final NCC abundance profile depends sensitively on the 
peak size r,jietai of the initial abundance profile. In runs D2 and 
D2-A, rjnetal is comparable to or smaller than the radius (~ 50-60 
kpc for A1795) within which AGN outbursts can efficiently mix 
metals. For these runs (D2 and D2-A), the central abundance 
peaks are effectively removed, result ing in flat NCC abundance 
profiles which are of ten observed (e.g. IDe Grandi &: Molendiir2001l : 
ILeccardi eFani^OTol l. 

the NCC to CC transformation, the increase of gas mass 
within ^ 60 kpc is less than one tenth of the gas mass 
within the original NCC metallicity peak. Consequently 
another powerful AGN outburst later on could not re- 
move this bro ad abu ndance peak. Furthermore, simula- 
tions by Poole et al.l (pOOS ) indicate that broad central 
abundance peaks can not be easily destroyed by mergers. 
Therefore, broad central abundance peaks, once formed, 
may persist for the remaining life of clusters. 

Some low-redshift CC clusters have narrow central 
abundance peaks (e.g., A2199). These systems may 
have stayed in the CC state for a relatively short dura- 
tion, otherwise a broad central abundance peak would 
have been established by AGN feedback. Interest- 
ingly, current X-ray observations seem to indicate that 
most low-redshift (z < 0.1; Dc Grandi ct al. 200 4 and 
interm ediatc-rcdshift (0.1 < z < 0.3; IBaldi et alTl2007l : 
ILeccardi fc Molendi 2008 ) CC clusters have quite broad 
central abundance peaks (~ 0.1-0.2r2oo)- This may sug- 
gest that most low-redshift NCC clusters without central 
abundance peaks were formed at higher redshift, when 
more CC clusters had smaller rmctai- Observers often 
emphasize m ean metallicity profiles averaged over many 
clust ers (e.g. . IDe Grandi et a l. 2004; Leccardi fc Molendil 
120081) ■ but our analysis suggests that it may be informa- 
tive to statistically study the distribution of rmctai in a 
large sample of CC clusters (and NCC clusters as well). 
An analysis of temporal evolution of the averaged rectal 
over cosmic time is also of great interest. The relative 
fraction of both types of NCC clusters and its redshift 
evolution may reveal the evolution history of the ICM 
thermodynamics and AGN activity. 

Alternativ e ly, r ecent cosmological s imula tions by 
iBurns et all (|2008l ) and iHenning et all (|2009f ) suggest 
that the origins of CC and NCC clusters depend on their 
different merger histories: CC clusters accreted mass 
slowly over time and grew enhanced CCs via hierarchical 
mergers, while NCC clusters experienced major mergers 
early in their evolution that destroyed embryonic CCs 
and produced conditions that prevent CC reformation. 
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Fig. 11. — Time evolution of radial profiles of emission- weighted 
spherically averaged gas entropy (top panel) and metallicity (bot- 
tom panel) in run Dl. Due to the powerful AGN outburst, the 
cluster reaches the NCC state at t ~ 0.2 - 0.3 Gyr. As the ICM 
cools at later times, it settles inward gradually. Ilowever, the re- 
sulting decrease in the size of the central abundance peak is not 
large (r^gtal ~ 130 kpc at t ~ 3 Gyr). 

However, as we discussed in Section 1, current cosmo- 
logical simulations do not have enough resolution to ac- 
curately model core structure and evolution. Further- 
more, ma ny NCC clu sters have a short cooling time 1 
Gyr; San derson et a l. 2006), and it is unclear why ra- 
diative cooling should not change their thermal state if 
they are formed early in their evolution history. Simi- 
larly, this scenario could not explain the origin of NCC 
clusters with flat metallicity profiles; central abundance 
peaks would be established by long-time metal enrich- 
ment from the BCGs if powerful AGN outbursts are not 
invoked to destroy them. 

4. DISCUSSION AND SUMMARY 

X-ray observations of galaxy clusters indicate a strik- 
ing bimodality in the properties of cluster cores, which 
separate clusters into two distinct classes: CC and NCC 
clusters. The origin of this bimodality remains un- 
clear. At the same time, a handful of extremely pow- 
erful, core-changing AGN outbursts have recently been 
dete cted in clus ters with a total en ergy ^ 10^^ — 10^^ 
erg (|Nulsen et al. 2005a,b; McNama ra et al.l l2005f). By 
conducting a suite of two-dimensional axisymmetric hy- 
drodynamical simulations, we inquire if such strong AGN 
outbursts can transform a CC cluster to the NCC state, 
and if so, how this transformation happens. Cluster 
abundance profiles have been well observed by X-ray ob- 
servations and contain important information on phys- 
ical processes in clusters. Consequently, we also follow 
the evolution of iron abundance profiles and study how 
they are affected by AGN outbursts during the CC to 
NCC transformation. 

We assume that AGN energy is injected into the ICM 
mainly in the form of CRs, and follow the co-evolution 
of these CRs with the cluster gas in our simulations. We 
find that when a large fraction of this energy is injected 
near the cluster center, strong AGN outbursts with en- 
ergy ~ 10^^ — 10^^ erg (see Table 1) can completely re- 
move cool cores, transforming the CC cluster (A1795 as 
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our fiducial cluster) to the NCC state. In view of the high 
central gas density in CC clusters, deposition of AGN 
energy within the cluster core may be common. The 
deposition of CRs produces shocks in the ICM, which 
propagate outward and heat the core gas to high en- 
tropies 100 kev cm^) very quickly (within ~ 10-50 
Myr in our main run Dl). The CRs also drive a strong 
gas outflow, producing a large, low-density cavity near 
the cluster center (with a radius of 80 — 100 kpc), 
which breaks up at time t ~ 200 Myr. During the break- 
up of the cavity, high-entropy thermal gas flows back to 
the cluster center and is efficiently mixed. The cluster 
relaxes to the NCC state with a central entropy ~ 100 
kev cm^ by time t ~ 200-300 Myr. 

Recent X-ray observations reveal several extremely 
powerful AGN outbursts with energy ~ 10^^ — 10^^ 
erg, inclu ding Hydra A (INulsen et all '2005a'), Hcr- 
cules A (iNulsen et al.l j2005bl) and MS0735 .6 -^742 1 
(|McNamara et al.l l2005t [McNamara et al.l I2OO90 . 
MS0735.6-h7421 hosts the most powerful AGN outburst 
currently known with an estimated AGN energy ~ 10^^ 
erg, but most of the energy seems to be deposited 
beyond 40 kpc (see IGu o & Mathews 2010 for a detailed 
modeling of this outburst), and the CC in this cluster, 
though significantly heated, is not completely removed. 
In contrary, the temperature profile in Hydra A is quite 
flat, suggesting that the CC in this cluster may have 
been removed by the observed powerful AGN outburst. 
Furthermore, the observed "Swiss-cheese-like" topology 
ass ociated with thre e pairs of X-ray cavities in Hydra 
A ([Wise et al.ll2007[ ) is also well reproduced in our run 
Dl-A, where CRs are injected into the ICM continuously 
by a (jet) source moving out from the central AGN. 
Hercules A hosts a young {t ~ 59 Myr) AGN outburst 
(~ 3 X 10^^ erg), which has signi ficantly disturbed the 
cluster core (INulsen et al.l l2005bf ) . While Hercules A 
is not fully relaxed, there is already so me evidence 
sugge sting that it may be a NCC cluster ([White et al.l 
Il997| ). Deep X-ray observations are needed to reveal its 



detailed thermal structure. 

During the transformation of CC to NCC systems, 
AGN outbursts efficiently mix metals in cluster central 
regions (within ^ 50 — 60 kpc ~ 0.03r2oo for A1795), 
where abundance peaks are usually seen in CC systems. 
The maximum abundance in the central peak thus de- 
creases, but the peak is not removed if its spatial size 
('"metal) is large. However, when rmctai is comparable 
to or smaller than the characteristic radius within which 
AGN outbursts can efficiently mix metals (as in the clus- 
ter A2199), the central peak is effectively removed and 
the resulting NCC abundance profile is relatively fiat. 
Our model naturally explains these two distinct types 
of NCC clusters seen in observations - those with and 
without central abundance peaks - both of which can 
develop following powerful AGN outbursts in CC clus- 
ters with different initial spatial abundance peak sizes. 

In contrast, the merger scenario for creating NCC clus- 
ters may not explain those NCC systems wi thout strong 
centra l abundance peaks. Simulations by iPoole etali 
(|2008f ) indicate that merging clusters which initially host 
central abundance peaks do not yield merger remnants 
with flat metallicity proflles (but also see Vazza et al] 
I2010D , though the dependence of their results on the size 
of the initial abundance peak remains to be explored. 
A statistical analysis of the peak sizes and their redshift 
evolution combined with a study of the possibly redshift- 
dependent relative fraction of both types of NCC clusters 
may shed interesting insights on the evolution history of 
the ICM thermal state and AGN activity in galaxy clus- 
ters. 
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